Implementation of a CNOT gate in two cold 
Rydberg atoms by the nonholonomic control 

technique 

E. Brion^, D. Comparat^ and G. HareP 
^ Laboratoire Aime Cotton, CNRS II, 
Batiment 505, 91405 Orsay Cedex, France. 
'^Department of Computing, University of Bradford, 
Bradford, West Yorkshire BD7 IDP, United Kingdom. 

February 1, 2008 

Abstract 

We present a demonstrative application of the nonholonomic control 
method to a real physical system composed of two cold Cesium atoms. 
In particular, we show how to implement a CNOT quantum gate in this 
system by means of a controlled Stark field. 

1 Introduction 

Quantum control is a very topical issue which bears relevance to many different 
fields of contemporary physics and chemistry, such as Molecular Dynamics in 
laser fields and Quantum Optics |21 El E El El ■ A few examples of control 
of the quantum state by adiabatic transport El , by unitary evolution 6 or by 
conditional measurements [7| El have been already proposed for the particular 
quantum system of atoms interacting with quantized electromagnetic field in a 
single-mode resonator. 

In parallel, a theoretical framework of quantum control has been built up. 
Four different types of problems have been identified in the literature [HI 
the control of pure state, the control of density matrix, the control of observable 
and, finally, the control of the evolution operator. Each of these problems 
can be formulated in the same way: the goal is to impose on the considered 
characteristics an arbitrarily chosen value. 

To achieve a control objective, one has to perturb the system, since its 
natural evolution usually results in too restrictive a dynamics. The control 
Hamiltonian H (t) thus comprises the unperturbed Hamiltonian Hq as well as 
M Hamiltonians of the form C„i {t) P„i , which stand for the interaction Hamil- 
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tonians of the system with M classical fields of controllable amplitudes Cm (t) , 



The functions {Cm (t)} play the role of the control parameters one has to adjust 
in order to achieve the desired control process. Any problem of control can thus 
be put in the following form: for a given physical system, specified by Hq and 

|Pm|, find the values of the control parameters {Cm (t)} which ensure that a 

specific characteristics of the system (quantum state, density matrix, observable, 
evolution operator) will take an arbitrarily prescribed value. 

Not all the objectives are always feasible. For example, the unitarity of the 
evolution operator for closed systems prevents the eigenvalues of the density 
matrix from changing through a Hamiltonian process of control. This kind of 
constraints is often referred to as kinematical constraints But there also 

exist dynamical constraints which stem from the algebraic properties of the 
Hamiltonians < Hq, Pm \- Indeed, the evolution operator 



where T denotes the chronological product, belongs to the Lie group obtained 
by the exponentiation of the Lie algebra generated by the operators \ iHo, iPm \ ■ 



The feasibility of a particular problem of control in a specific physical situation, 
defined by the Hamiltonians |ifo,-Pm|j is clearly related to the properties of 
this algebra. For example, if one wants to completely control the evolution oper- 
ator of a quantum system, i.e. to be able to give the operator U any prescribed 
value, one must perturb the system in such a way that the operators \ iHQ, iPm \ 



generate the whole Lie algebra u {N) which provides, through exponentiation, 
the whole Lie group U [N] ^1 ^] (this prescription is called the Bracket Gen- 
eration Condition). Necessary mathematical conditions also exist for the other 
types of control problems which can be found in the literature 10 : these con- 
ditions are obviously weaker than the previous one, since the controllability of 
the evolution operator automatically implies all the other ones. 

The feasibility of a control problem can thus be decided through mathe- 
matical criteria established in the context of the Lie group theory. But the 
explicit values of the control parameters achieving the desired control objective 
still remain to be found. In other words, once the existence of a solution has 
been proved, one still has to find it explicitly. To achieve this goal, different 
methods, such as optimal control |14[ll5l l^. have been proposed, most of which 
rely either on a known or intuitively guessed particular solution which can be 
further optimized with respect to a given cost functional, through variational 
schemes [5]. A purely algebraic approach based on the decomposition of 
the arbitrary desired evolution on the Lie group, is also possible, but rapidly 
leads to intractable computations as the dimension of the state space increases. 
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In the context of the control of the evolution operator, a constructive method 
called nonholonomic control |17| was proposed, in the same spirit as in |18| : 
this method is fundamentally algebraic but also uses optimization steps. The 
physical idea is to alternately apply two distinct well chosen perturbations Pa 
and Pb {i.e. two perturbations which check the Bracket Generation Condition). 
The timings of the interaction pulses play the role of control parameters and 
are determined by solving the " inverse Floquet problem" . The convergence of 
the algorithm results from an unsuspected simplification which emerges from 
the Random Matrix Theory. Indeed, it relies on the algebraic properties of the 
]\j-th j-QQ^g Qf ^jjg identity matrix, the spectra of which resemble those of random 
unitary matrices which are ruled by the Dyson distribution law. This method, 
combined with a generalization of the Quantum Zeno effect, has also led to 
coherence protection schemes [T^l^ . 

The nonholonomic technique is completely universal and can thus be applied 
to any physical system. In the last few years, cold atomic Rydberg states have 
appeared as particularly relevant in the context of quantum information, and 
have been widely investigated both theoretically [21 [22] and experimentally 
[211 [Ml- In this paper we propose to consider a real system composed of two 
cold Cesium (Cs) atoms interacting via dipolar forces and prepared in Rydberg 
states. Under some physical assumptions, the Hilbert space of the system is 
restricted to four states, as for a two qubit system. We then propose a control 
experiment employing a pulsed electric field which allows us to impose a CNOT 
gate to the system through nonholonomic control. Even though we have tried 
to propose a realistic experimental setting, the simplifications we have made 
result in serious limitations of the scheme we present. Nevertheless, it shows 
how the nonholonomic method can work on a real physical system and suggests 
that this technique can be an effective way to solve real-life problems of control. 

The paper is organized as follows. In Sec. II, we recall briefly the main fea- 
tures of the nonholonomic control technique. In Sec. Ill, we describe our physical 
application in detail - the system, the control Hamiltonians, the calculated con- 
trol parameters - and discuss its validity and limitations. 

2 Control of the evolution through nonholonomic 
control 

Let us consider an A^-dimensional quantum system of unperturbed Hamiltonian 
Hq. Our goal is to control its evolution operator U, i.e. to be able to achieve any 
arbitrary evolution Uarbitrary G U{N). To this end, we alternately apply two 
physical perturbations, of Hamiltonians Pa and Pb, during iV^ time intervals 
(pulses), the timings of which are denoted by {rfc =tk — ife-i}fc=i jv2 (io = 
and tjqi = T correspond to the beginning and the end of the control sequence, 
respectively). The control Hamiltonian takes the following pulsed shape (cf Fig. 

m 

Hit) ^Ho + Ca it) Pa + C'b (t) Pb 
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Figure 1: Pulsed shape of the control Hamiltonian. 



with Ca [t) = 1, Cb [t) ^ and H{t) = Hq + Pa ^ Ha for t e [tafe, fsfc+i] and 
Ca (t) = 0, Cb (t) = 1, and H{t) ^ + Pb = Hb for t G [t2k-ut2k\ , where 
fc = 1, . . . , iV^. The total evolution operator is therefore 

& ({n, . . . , rjv2^i, rjv2}) = e-^-^^^«^ • e-^^""«^-i . . . e"*^"^! , 

where we have implicitly assumed that N is even. 

Our control problem can thus be translated into the following form: given 
an arbitrary unitary operator U arbitrary G U{N), we want to find a time vector 
T = (ri, . . . , Tjva) with non- negative entries, such that 



As we said previously, for a solution to exist the operators |i-ffa, iHb | must 
generate the whole Lie algebra u (N) . This prescription can be checked directly 
as long as the dimension N is not too big: one simply computes the commutators 
of all orders of iHa and iHi, and stops as soon as they generate u{N). But 
when N becomes large, direct computation is intractable. In that case, one 
can simply check the following sufficient conditions, suggested by V.G. Kac 
|17| . according to which the system becomes nonholonomic, that is completely 
controllable, when the representative matrix of Hb in the eigenbasis of Ha has 
no zero elements and vice versa, and also the eigenvalues and their pairwise 
differences are distinct for both matrices. 

Once the previous criterion is checked, one has to compute the time vector 

solution of Eq.l^. The method consists in determining the time vector 
such that U {'t^^^'') ~ I, then iteratively approaching the time vector V through 
a Newton-like technique. 

The straightforward way to compute T**^"^ would be to minimize the func- 
tional F {t ) — ||[/ ( T ) — with respect to r . However, F presents many 



U (V) — Uarbit 



.rary • 



(1) 
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local minima which make its optimization uneasy. Yet, there exists an alterna- 
tive method based on the algebraic properties of the TV*'' roots of the identity 
matrix. 

The idea is to look for N parameters {Tfc}fe=i..jv such that 



(2) 



where /« is a non-degenerate N*^ root of the identity matrix, i.e. a matrix of 
the form 

" 1 ••• 
e*T^ •■• 



In = M- 







•M, 



where M is a unitary matrix. To compute the Tfc's, we use the following al- 



gebraic property: if Pfyi^) = X]iLo%'^"' denotes the characteristic polynomial 



of a unitary matrix U, then ^ 



N I I 



> 2 and the equality is achieved if and 



only if U is an N*^ root of the identity matrix, up to a global phase factor. 
To obtain the T^'s, one thus computes the characteristic polynomial 



N 



Pp(A)^5]a,({Tfe},..i...^)A^" 



of the matrix product 



. . e 



and minimizes the function 



N 



.N 



to 2 with respect to the T^'s. This minimization turns to be quite easy, due 
to the fact that a generic unitary matrix is very close to an iV*'' root of the 
identity. In fact, numerical work shows that in about 30% cases of randomly 
chosen timings {Tk } the standard steepest descent algorithm immediately finds 
the global minimum Fn = 2. This fact has roots in the Random Matrix Theory 
pS) . Indeed, according to Dyson's law, the eigenvalues of random unitary ma- 
trices tend to repel each other, and are thus very likely to be almost regularly 
distributed on the unit circle, as those of an iV'^ root of the identity, as shown in 
Fig[21 In other words, in the space of iV x TV unitary matrices, the /« matrices 
are present in abundance, and can be reached from a randomly chosen point by 
small variation of the timings. 

Finally, we define the time vector l^^") corresponding to the identity matrix 
by simple repetition of {Tk} 



\+U-l)N 



T, for 



(3) 
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Figure 2: Spectrum of a random unitary matrix (white circles) resulting from 
the repulsion of the eigenvalues on a unit circle is shown vs the eigenvectors 
of -/V*'' root of the identity matrix (black circles) multiplied by an unimportant 
phase factor. 

and checks that indeed 



up to an irrelevant global phase factor. 

We now have to iteratively determine the time vector V from ~t^^^K Let us 
first consider the case of a target evolution close to the identity: in that case, 
Uarbitrary Can bc Written in the form 

Uarbitrary = C/e = exp(-lHe), (4) 

where Ti is an NxN bounded (| |7i 1 1 < 1) dimensionless Hermitian Hamiltonian, 
and e > a small parameter. We then calculate the variations St^, which are 
determined to first order in e by the linear equations 




Once (5 V = {Sti , . . . , St[^2 ) has been calculated through standard techniques of 
linear algebra, we replace by V^"' + and repeat the same operation 

until we obtain t which checks U ( t ) = Uarbitrary at the desired accuracy. 



N ^ 



If the evolution UarUtrary — IS not close to the identity, that is, if e 
is not small, one has to divide the work into elementary paths on which the 
previous method converges. To this end, we consider an integer n > 2 such that 

(^Uarbitrary^ = [/jL is attainable from / through our iterative algorithm, and 
determine in this way the associated time vector ^(i) which checks 

1 

Taking {UarUtrary^ = ' ^ as our new target, we repeat the same algorithm 
to compute t such that 

and so on. We progress in this way as long as the algorithm converges: in 
general, it stops at a value n* > 1, for which Eq.© has no solution. Then, we 
keep the time vector and simply repeat the same control sequence n* 

times to achieve the desired evolution 

^(^(-.))---^(^(^)) - 

— U arbitrary- 

To conclude this section, let us point out that another equivalent form of 
the nonholonomic control method exists, in which the durations of the pulses 
are fixed, whereas the strengths of the perturbations play the role of control 
parameters. The equations in that case are very similar to those we have just 
dealt with, and the same method applies with almost no change. For more 
details, see [T7] . 



arbitrary 



3 Application: implementation of a CNOT gate 
in a two cold Cs atom system 

In this section, we present the application of our control technique to a system 
of cold Cs atoms. Frozen Rydberg gases of interacting Cs atoms have been 
investigated in |26| and have revealed the existence of new phenomena typical of 
low temperatures, such as the modification of resonance profiles, the explanation 
of which requires the framework of a TV body theory. The system we have chosen 
to consider in this paper is greatly inspired by the experimental situation studied 
in it consists in two Cs atoms in Rydberg states, denoted by {A) and {B), 

of dipole momenta (Ia and ds , respectively, linked by the fixed vector R = R~n 
which is determined by its norm, taken equal to R = (of the same order 
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as the distance between two close neighbour atoms in [221 )j and its direction 
it, defined by polar angles 9 and ip (cf Figl^J. These two atoms are coupled by 
dipole-dipole interaction 



1 

Vdd 



3{dA-n]{dB-n 



and are subject to a Stark field £'5 — Es e z, the z axis corresponding to 
the quantization axis for the total angular momentum. The total Hamiltonian 
of the system is thus composed of the unperturbed part, Hq = Vdd, and the 

controllable perturbation Vs = — d ■ Es- 

At this stage, we shall make two remarks. Firstly, for the dipole-dipole 
approximation of the interaction energy to hold, i.e. for higher order terms to 
be actually negligible, the distance R between atoms must be much greater than 
the sum of the radii of the atoms, which is not strictly the case here: indeed, 
in the atomic Rydberg states we shall consider {n = 23,24), the two atoms 
have almost the same radius, approximately equal to 3 ■ 10~^m, whence the 
prescription R ^ ta+tb is not rigorously checked. To circumvent this difficulty, 
one could be tempted to increase the value of the interatomic distance R, but this 
would result in a decrease of the typical value of the dipole-dipole interaction 
which would then become much smaller than the typical value of the Stark 
interaction energy: this has to be avoided for our purpose of control, and the 
balance between the unperturbed Hamiltonian Hq = Vdd and the perturbation 

Vs = — d ■ E s must be preserved. One might then suggest to decrease the Stark 
field in order to make the Stark interaction term decrease and follow the dipole- 
dipole interaction, but then the range of the required values for the Stark field 
would not be realistic. We choose instead to keep the value of i? = 2 • 10~^m and 
to consider the dipole-dipole term only, bearing in mind that a more rigorous 
approach should take higher order effects into account. Secondly, contrary to 
the experimental situation described in |2fij . we do not deal with a sample of 
A'' interacting Cs atoms but only with a two atom system. Nevertheless, it has 
been demonstrated that this situation could be experimentally achieved |27|. 

Let us now describe the control operation we want to achieve. Initially in 
zero field, the system is prepared in an arbitrary superposition J2^=o N) 
following four states 

|0) = |24si/2,m, = l/2;23si/2,mj = 1/2) 

|1) = |23p3/2,"ij = 3/2;23p3/2,mj = 3/2) 

|2) = |23p3/2,m, =3/2;23p3/2,mj = 1/2) 

|3) = |23p3/2,m, = l/2;23p3/2,mj = 1/2) 

which formally stands for two qubits of information to be processed. Our goal 
is to impose the CNOT gate, or, in other words, whichever the initial state 

1(^(0)) = CO |0) + ci|l)+ C2|2)+C3|3) 
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Figure 3: Polar angles 9, ip. 



is, we want to impose the particular evolution CNOT, yielding the final state 

CNOT Iv? (0)) = Co |0) + ci |1) + C3 |2) + C2 |3) . 

Let us underline that the choice of this specific gate is purely arbitrary: we 
could have chosen any other unitary evolution of the system. Nevertheless, the 
CNOT gate is particularly important since it enters into the composition of 
universal sets of quantum gates, as demonstrated by D.P. DiVincenzo in his 
well-known paper JSH: ■ 

The control sequence which allows us to achieve this objective roughly con- 
sists in alternately and diabatically (i.e. abruptly) applying the two different 
values Ea = 87. 42^/ cm and = 84. 85^/ cm of the Stark field to the sys- 
tem, corresponding to two different values Ha and Hi, of the total Hamiltonian, 
during pulses the durations of which will be determined through the algorithm 
described in the previous section. From now on, we shall make the two follow- 
ing physical assumptions on the system: (i) the states |23si/2,mj = 1/2) and 
|24si/2,"^j = 1/2) are not mixed with the Stark multiplicities n = 19, 20 (this 
assumption is motivated by the position of the avoided crossings, which are far 
from the values Ea and E^ of the applied field); [ii) the states |23p3/2, mj = 1/2) 
and |23p3/2,?7ij =3/2) remain unaltered (we neglect their mixing with 22d 
states, which is correct up to 10%), while their energies decrease linearly with 
the amplitude of the applied field. According to these simplifications, the 
spectrum of our system in a static electric field can be represented as shown 
on Fig^ the energies of the states |1),|2) and |3) vary linearly with the 
appfied field (with the same slope 7 = —283.044 atomic units), while the 
energy of the state |0) remains constant; moreover, for the resonance fields 
El = 88.8V/cm = 1.73- lO-^a.u.,E2 = 84AV/cm = 1.64- 10-*^a.u., = 
8Q.5F/cm = 1.57- IQ-^a.u., we have E = E{\0)), E {\2)) = E {\0)) and 
E (|3)) = E (10)), respectively. In the basis {[0) , |1) , |2) , |3)}, the total Hamil- 
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Figure 4: Simplified Stark diagram of the system considered: with (dashed hne) 
and without Vdd (sohd hne). 



tonian takes thus the foUowing expression 
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which wiU alternately take the two distinct values Ha for E ^ Ea and Hi, for 



In summary, we deal with a system whose Hilbert space is restricted to 
iV = 4 states, to which we want to apply the CNOT gate. To this end, we 
propose to alternately apply two Hamiltonians Ha and Hb during pulses whose 
timings {t^} are to be determined by the method we presented in the previous 
section. What we have to do first is to find the iV = 4 timings {Tk=i^...,4} 
which meet Eq.Q by minimizing J2'j=oWji{'Ek=i,...,4})f- Using Eq.® we 
then build the = 16-dimensional time vector from the Tfc's, which 

achieves the evolution /. Then we apply the iterative algorithm described in 
the previous section. As the CNOT gate is far from the identity matrix, we 
have to divide the work: we take CNOT~ as our target evolution, where n 
is an integer greater than 1 for which our algorithm converges, providing the 
time vector V" ; then we take the new target CNOT'^ where n' < n and 
run our algorithm again, yielding the vector r ^ , etc. as long as we obtain 
convergence. The smallest value of n we obtained is n* = 8, associated with the 
time vector r ~ which achieves the evolution CNOT~ : the desired evolution 
CNOT is obtained by repeating n* = 8 times the same elementary control 
sequence, defined by We thus see that the time vector '^c:^OT 

achieves CNOT is n* x N^ — 8 x 16 — 128-dimensional and can be built by 
repeating 8 times the vector V^. 



E = Eb. 
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6, 77ns 



Figure 5: Elementary iV^ ~ 16 pulse control sequence obtained for R = 2 ■ 
lO^^m, = 7r/15, and = tt/6. The total duration is 6,77ns and the values 
of the different interaction times are ti{ns) = {0.743378, 0.553823, 0.277301, 
0.133699, 0.800748, 0.423586, 0.427981, 0.0427037, 0.71635, 0.458626, 0.403841, 
0.13241, 0.682599, 0.54579, 0.349389, 0.0809227}. 



Fig[Sl presents the numerical results we obtained for 9 = 7r/15 and (j) — 7r/6 
(cf Fig|3J) which shows the switchings of the Stark field on an elementary control 
sequence, whose duration is 6.77ns. In our calculations, we tried to remain in a 
realistic range for the different parameters of our system. For example, the total 
control duration, which is of the order of 8 x 6.77ns ~ 0.06/is, is much smaller 
than the lifetime of the Rydberg states considered, which is approximately 10/is. 
Yet, serious problems and limitations arise. 

Firstly, the required switching time of the Stark field is much less than 
Ins (some timings are a few lOOps), which is experimentally very difficult to 
achieve: this will unavoidably threaten the reliability of the control. To address 
this problem, one may consider replacing the static control fields by pulsed 
lasers, as in |19| . which would probably allow rapid switching times and would 
certainly be more tractable experimentally. Secondly, the four state system we 
have considered here is a severe idealization: the couplings between the states 

|24si/2,mj = ±l/2;23si/2,mj =±1/2) 
|23p3/2,TOj = ±3/2; 23p3/2, mj = ±3/2) 
\23p3f2,mj = ±3/2;23p3/2,mj = ±1/2) 
\23p3/2,m, = ±l/2;23p3/2,"^j - ±1/2) 

cannot be ignored. In addition, the influence on the states 23s and 24s of the 
multiplicities n = 19, 20 has been completely neglected, as well as the mixing of 
the states 23p3/2 with the states 22d. One can solve these problems by increasing 
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the state space, i.e. by taking all the states which are actually coupled by 
the Stark field and the dipole-dipole interaction into account: calculating the 
control time vector becomes much longer, as the system considered is larger, 
but, fundamentally, the structure of the problem remains the same. Finally, if 
we do not work with two atoms but rather with a large sample, it might be 
experimentally very difficult to fix precise values to R,6 and ip: this results 
again in a loss of reliability of the control. A possible solution to this problem, 
though not perfect, would be to put the atoms in an optical lattice, which would 
allow one to control more precisely their spatial arrangement. Another method 
would be to perform a first control sequence, the goal of which would be to 
distinguish between "good" and "bad" pairs of atoms: for instance, starting 
from the state |3), the "good" pairs {i.e. the pairs with the required vector 
R) will undergo the CNOT gate and will thus end in the state |4), whereas 
the other pairs will end in a superposition of all states and could therefore be 
experimentally distinguished and destroyed. 

To conclude this section, we want to emphasize that the limitations disc;usscd 
above do not remove the pedagogical and demonstrative value of the application 
presented. The example considered here shows the operability of the nonholo- 
nomic control method and suggests that it can be actually employed to achieve 
real objectives of control. 

4 Conclusion 

In this paper, we first recalled the nonholonomic control technique, which al- 
lows one to control the evolution operator of generic quantum systems which 
meet the bracket generation condition: after putting it in the general frame- 
work of quantum control we briefly exposed its main algorithmic features and 
underlined the fundamental reasons for its convergence. Then we suggested a 
demonstrative application of this scheme to a system of two cold Cs atoms, in- 
spired by experimental studies on cold Rydberg gases: we showed that through 
alternately applying two different values of a Stark field during 128 pulses, the 
timings of which range from AOps to SOOps, one can impose the CNOT gate to 
two qubits of information stored in four specific states of the system. Finally 
we discussed the physical validity and the limitations of our application. 

The authors thank V.M. Akulin and P. Pillet (Laboratoirc Aime Cotton, 
Orsay, Prance) for stimulating and fruitful discussions. The support of EU 
(QUACS RTN) is kindly acknowledged. 
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